Thermal Conductivity and Specific Heat of the Spin-Ice Compound Dy 2 Ti 2 07: 
Experimental Evidence for Monopole Heat Transport 
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Elementary excitations in the spin-ice compound Dy 2 Ti207 can be described as magnetic 
monopoles propagating independently within the pyrochlore lattice formed by magnetic Dy ions. 
We studied the magnetic-field dependence of the thermal conductivity for B\ \ [001] and ob- 

serve clear evidence for magnetic heat transport originating from the monopole excitations. The 
magnetic contribution fv mag is strongly field-dependent and correlates with the magnetization M(B). 
The diffusion coefficient obtained from the ratio of K mag and the magnetic specific heat is strongly 
enhanced below 1 K indicating a high mobility of the monopole excitations in the spin-ice state. 

PACS numbers: 66.70.-f, 75.40.Gb, 75.47.-m 
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The recent prediction of magnetic monopoles in the 
spin-ice compounds has attracted a lot of interest [D- 
llO] . Spin ice is a geometrically frustrated spin sys- 
tem, which is realized in Dy 2 Ti 2 07 by a sublattice of 
corner-sharing Dy 3+ tetrahedra. Due to a strong Ising 
anisotropy, the magnetic moments of Dy 3+ point either 
in or out of a tetrahedron. The magnetic dipole en- 
ergy is minimized, when two spins point in and two out 
of a tetrahedron (2in-2out), what is realized by 6 out 
of 2 4 = 16 possible configurations for a single tetrahe- 
dron. In this respect, the spin orientation of Dy 3+ cor- 
responds to the hydrogen displacement in water ice 11 1 



and the ground state is highly degenerate with a residual 
entropy S = Nk B /2 ln(3/2) for T->0K [TM^. Ex- 
cited states can be created by flipping one spin, resulting 
in two adjacent tetrahedra with configurations 3in-lout 
and lin-3out, respectively. In zero magnetic field, such 
a dipole excitation can fractionalize into two individual 
excitations, a monopole (3in-lout) and an anti-monopole 
(lin-3out), which can propagate independently. This can 
be visualized by flipping, e.g., another in-pointing spin 
of the 3in-lout tetrahedron such that it relaxes back to 
(another) 2in-2out ground state configuration, while the 
3in-lout state has moved to a neighboring tetrahedron. 

The model of magnetic monopoles has been widely 
used to describe many experimental observations of 
Dy 2 Ti907fl1-[loj|. Nevertheless, there are basic properties 
of the spin-ice materials, which are far from being under- 
stood. For example, the specific heat c p of Dy 2 Ti2 07 
has a pronounced maximum around 1.2 K resulting 
from the magnetic excitations and the corresponding 
entropy is close to the expected So = 1.68 JK _1 mol^y 
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lalHl- Below - 600mK, however, c p {T) data 



published by several groups differ by almost an order 
of magnitude, see figure 

ED H E3-G3- Very recently it 
became clear that in this low-temperature regime the 
magnetic subsystem of Dy 2 Ti207 enters a region of very 
slow dynamics with relaxation processes that may ex- 
tend over extremely long time scales [l(J E3] ■ Another 
open issue is the dynamics of the magnetic monopoles. 



The possible observation of a monopole current in an 
external mag netic field is currently under strong de- 
bate [1, 0, [l9| • In this context, it is also unclear whether 
the magnetic monopoles contribute to the energy trans- 
port and how the monopoles interact with each other and 
with the phonon excitations. A suitable probe to study 
these issues are measurements of the thermal conductiv- 
ity k, which are the main subject of this report. Only 
one experimental study about the heat conductivity of 
Dy 2 Ti 2 C>7 has been published so far 17] . This previ- 



ous work focuses on the anomalously enhanced relaxation 
times of Dy 2 Ti 2 07 in the low-temperature region, which 
are analyzed by assuming that the thermal conductivity 
is of purely phononic origin. In this report, we present 
a detailed study of the magnetic-field dependent thermal 
conductivity k(B), which yields clear evidence that up 
to almost 50% of the low-temperature heat transport of 
Dy 2 Ti 2 0*7 is of magnetic origin and our analysis suggests 
that this is a consequence of the high mobility of mag- 
netic monopoles in zero field. 

For the actual study, large single crystals of Dy 2 Ti2C>7 
and Y 2 Ti 2 07 were grown by the floating-zone technique. 
The crystals were oriented in a Laue camera and samples 
of typical sizes of several mm 3 were cut. The thermal 
conductivity was measured by the standard steady-state 
method from ~ 0.3 to 300 K. The temperature gradi- 
ent was produced by a heater at one end of the sample 
and measured by two matched Ru02 thermometers for 
T < 15 K in magnetic fields up to 1 T. In a separate run, 
the measurements were extended up to 300 K by using 
AuFe(0.07%)-Chromel thermocouples to detect the tem- 
perature gradient. The magnetic field was applied along 
[001], perpendicular to the heat current along the [110] 
direction of a sample of 1 x 1 x 3 mm 3 with the long edge 
parallel to [110]. For this sample geometry, demagneti- 
zation corrections have to be taken into account. The 
magnetization, needed for the correction, has been mea- 
sured with a home-built Faraday magnetometer on a thin 
sample of 0.3 x 2 x 4 mm 3 with the long axis || [001] to 
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FIG. 1. (color online)- Comparison of our specific heat data 
with literature data [21. Ilq. 1 1 71] - The low-temperature plateau 
of our data (open symbols) probably arises from nuclear con- 
tributions of 161 Dy and 163 Dy. Inset: Schematic illustration 
of the method to measure the specific heat. 



minimize demagnetization effects. 

The specific heat was measured in a home-built 
calorimeter using a quasi-adiabatic heat-pulse method. 
A sample of 21 mg was fixed to the platform by a small 
amount of grease. The addenda was measured in a 
separate run and was subtracted. It has been shown 
recently 



17] . that below 600 mK the temperature re- 



laxation of Dy 2 Ti207 contains multiple time constants, 
which can be understood by Dy 2 Ti207 consisting of sub- 
systems weakly coupled to each other and to the plat- 
form. Standard methods to measure c p do not account 
for such sample-internal dynamics, and thus provide too 
small values [13] ■ Our approach to measure c p below 
600 mK, slightly differs from that of [17]. As the sam- 
ple is at low temperature and in high vacuum we can 
safely assume that the subsystems C, are not directly 
linked to the bath, but only to each other and to the 
platform (inset of figure Q}, whose temperature Tp is 
measured as a function of time. Starting from equilib- 
rium (T samp i = Tp = Tbath), a constant heating power 
P is applied to the platform until saturation is reached 
(T samp i c ~ T P = T bath + AT [20]). The heat AQ stored 
in the sample is then calculated from the total dissipated 
heat by subtracting the numerically obtained heat flown 
from the platform via K to the bath. Above 600 mK, 
this method yields heat capacities consistent with those 
obtained by conventional single-relaxation methods. At 
lower temperatures, however, this technique yields en- 
hanced Cp(T) data, similar to those of |17j, but with- 
out any a-priori assumptions about the subsystems or 
their respective couplings. Figure [1] compares our re- 
sults with literature data [1, 16, 17[, which all agree with 
each other within 10% above 600 mK. For lower temper- 
atures, however, the standard techniques [H, [lj| result in 
significantly lower c p values than those obtained by the 
methods which explicitly account for the glassy behav- 
ior of different subsystems in Dy 2 Ti2 07. In comparison 
to fl7j j . our data show an additional shoulder in c p (T) 
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FIG. 2. (color online). Thermal conductivity of Dy 2 Ti207 for 
zero field and for 0.5 T and the zero-field thermal conductivity 
of the non-magnetic Y2Ti207. Inset: Same data on linear 
scales up to 250 K. 



around 500 mK and a tendency towards saturation be- 
low 350 mK. We suspect that nuclear contributions to 
c p resulting from the isotopes 161 Dy and 163 Dy [2l[ may, 
at least partly, cause the anomalous low-temperature be- 
havior of our data, which is not observed in 17J, where 



a b Dy-enriched sample was studied. In order to clar- 
ify the origin of these deviations, further studies on dif- 
ferent samples are necessary. For the following discus- 
sion, this uncertainty is of minor importance, because (i) 
the thermal conductivity k(B) has only been measured 
above 350 mK and (ii) the deviation between our c p (T) 
and the data of [ItJ around 500 mK will affect only a 
single point of the diffusion coefficient discussed below. 
Nevertheless, we emphasize that due to the glassy low- 
temperature behavior of Dy 2 Ti2 07 below 600 mK any 
quantitative comparison of experimental data with theo- 
retical models should be treated with some caution. 

The thermal conductivity k(T) of Dy 2 Ti2C>7 was mea- 
sured at zero and in an external magnetic field of 1 T 
(figure [2|) , but as mentioned above the external field was 
corrected by a considerable demagnetization field of up 
to 0.5 T. Above 4K, there is almost no field-dependence 
of k(T), while k is significantly suppressed by a magnetic 
field for lower temperatures. In order to obtain informa- 
tion about the phononic contribution of k(T), we also 
studied the iso-structural, but non-magnetic Y2Ti2C>7. 
As can be seen in the inset of figured! k(T) of Dy 2 Ti207 
and of Y2T12O7 are very similar above 100 K, but at 
lower temperature k(T) of Y2Ti2C>7 is significantly higher 
than that of Dy 2 Ti2C>7. Probably, this difference is re- 
lated to an anomaly observed at 110 K in Raman spec- 
troscopy data that indicates a structural instability in 
Dy 2 Ti 2 Cv which is absent in the non-magnetic homo- 
logue Lu 2 Ti 2 7 [H, El- In addition, « of Dy 2 Ti 2 7 
can be reduced by phonon scattering via crystal-field ex- 
citations of the partly filled 4/ shell of Dy. Consider- 
ing the low-temperature range, it turns out that n(T) 
of Y 2 Ti 2 7 follows a power-law k(T) oc T 2A below 3 K 
and a similar behavior (k(T) oc T 2 2 ) is present for the 
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k(T) data of Dy 2 Ti207 in a magnetic field of 0.5 T. In 
contrast, the zero-field k(T) of Dy 2 Ti 2 C>7 shows a clear 
shoulder around 1 K. This qualitative difference sug- 
gests the existence of an additional magnetic contribu- 
tion K ma g, which appears in the zero-field data on top of 
the phononic background, that is 
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where K p h is essentially represented by the k(T) data 
measured in 0.5 T. This conclusion is confirmed by mea- 
surements of the field-dependent k(B) at constant tem- 
peratures. Again, demagnetization has been taken into 
account to rescale the magnetic field. Figure [3] displays 
the relative change k/kq for different constant tempera- 
tures. Below 2 K, we find a step-like decrease of k(B) 
around 0.2 T, which systematically sharpens on decreas- 
ing temperature. As is shown exemplarily for 800 mK in 
figure OJi, this step is followed by a weak, essentially lin- 
ear decrease towards higher fields. The relative reduction 
k/kq has a maximum around 600 mK and vanishes prac- 
tically above 4 K. The decrease of k{B) correlates with 
the increase of M(B), shown exemplarily for 700 mK and 
500 mK in figure [3p and [3ji, respectively. The step-like 
change of k/kq exhibits a clear hysteresis below 700 mK 
(figure [3j>e) , with different critical fields depending on 
the field-sweep direction. 

Finally, an additional feature appears in k(B) below 
500 mK (figure HJs) : Measuring an initial curve (1) after 
zero-field cooling and subsequently decreasing the mag- 
netic field back to zero (curve 2), fc(0T) only recovers 
about 90% of its original zero-field value. Repeating 
the field sweeps from this new starting point results in 
k(B) curves (3) and (4) with coinciding endpoints, where 
curves (2) and (4) perfectly match each other. The re- 
duced zero-field values k(0 T) slowly relax back to the re- 
spective zero-field-cooled values, as is shown in Figure Ef 
for 400 mK. In order to describe this relaxation process 
at least two relaxation times are needed, that is 

«(t) = ao + oi(l - e-*/ Tl ) + 02(1 - e-*/ T2 ) . (2) 

The fit of figure [3f yields large relaxation times 
T\ ~ 8 min and r 2 ~ 100 min. As every data point of k(B) 
in figure[3K-e requires several minutes of temperature sta- 
bilization, the different zero-field values (figure Bfe) orig- 
inate from the extremely slow relaxation r 2 [24]. This 
very slow relaxation only occurs after a field sweep be- 
low 500 mK but not after cooling at zero field. Such a 
slow glass-like low-temperature behavior of Dy^Tio O? is 

in agreement with various recent reports slums. 

A straightforward qualitative interpretation of the ob- 
served magnetic-field dependent n{B) is obtained by as- 
suming that the magnetic contribution K mag results from 
a heat transport by magnetic monopoles. Due to the de- 
generate zero-field ground state (2in-2out), a monopole 
excitation (3in-lout or lin-3out) can easily propagate by 
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FIG. 3. (color online). Field dependence of k(B)/k(0T) for 
B|| [001]. All curves were measured after zero-field cooling. 
The fit of the linear decrease above ~ 0.3 T is shown exem- 
plarily for 800 mK. Panels (b) and (d) also contain the mag- 
netization data M(B) (A), and (f) displays the time depen- 
dent relaxation k(0T, t) back to the initial zero-field value for 
400 mK. 



single spin flips. When a monopole excitation passes 
through a 2in-2out tetrahedron, two subsequent single 
spin flips are needed, which change the initial 2in-2out 
configuration of this tetrahedron to another 2in-2out 
state. As in zero magnetic field all 2in-2out configu- 
rations are degenerate, the magnetic monopoles have a 
large mobility. This drastically changes in a magnetic 
field along [001], which lifts the degeneracy of the differ- 
ent 2in-2out configurations and thus the monopole mo- 
bility is suppressed. In addition, the correlation between 
increasing magnetization and decreasing heat conductiv- 
ity is explained. The magnetic field determines, which 
2in-2out state is the field-induced ground state, and for 
not too small fields (> 0.1 T), the magnetization is pro- 
portional to the population of this particular 2in-2out 
state, whereas the monopole mobility systematically de- 
creases the more this 2in-2out state is populated. With 
increasing £?||[001], the monopole excitation energy also 
increases weakly, i.e. the monopole density decreases, 
but we expect this effect to be of minor importance, be- 
cause the decreasing monopole mobility due to the field- 
induced lifting of the ground-state degeneracy will be the 
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FIG. 4. (color online). Magnetic contribution fi: mag of the 
zero-field heat transport (inset; the line is to guide the eye) 
and the diffusion coefficient -D mag = K ma g/c ma g, together with 
the estimated mean free path I, 

dominating effect. For magnetic fields above the step-like 
decrease of k(B), the heat transport is purely phononic 
and in order to extract the zero-field phononic contri- 
bution, we extrapolate the weak linear decrease of n(B) 
back to T (figure . The extrapolation uncertainties 
determine the error bars of K mag . The resulting zero- field 
magnetic contribution n mag has a pronounced maximum 
at ~ 1.2 K (inset of figure 2]). The maximum of K mag is 
located close to that of c mag , but on further increasing 
temperature K mag rapidly decreases and practically van- 
ishes above ~ 4K. The main panel of figure [4] displays 
the corresponding diffusion coefficient, 



Cmag 

Above 1 K, -D m ag is almost temperature independent but 
strongly increases towards lower temperatures. This fits 
to the qualitative expectation of a high monopole mobil- 
ity by single spin flips on the degenerate ground state. At 
low temperatures, the monopole excitations are highly 
dilute resulting in a large mean free path and hence a 
large diffusion coefficient. With increasing temperature, 
the number of monopoles increases and thus the diffu- 
sion coefficient is expected to decrease. In the simpli- 
fied model of single spin flips, this follows from the fact 
that a monopole can pass an already excited tetrahedron, 
i.e. another (anti-)monopole, only via simultaneous spin 
flips. Within kinetic gas theory, the diffusion coefficient 
is related via D = v£/3 with the mean velocity v and 
the mean free path I of the particles, but this relation 
cannot be directly applied to the actual monopole gas 
in Dy 2 Ti2C>7 for various reasons. First of all, the num- 
ber of monopoles is not conserved and, moreover, neither 
the (average) velocity of monopoles nor their interaction 
with each other or with the phonon excitations is well 
understood up to now. Quite recently, an expression for 
the monopole mobility in a (magnetic or electric) field 
has been derived 0] , but it is unclear whether this result 
can be related to our experimental -D mag arising from a 



finite temperature gradient in zero field. Despite these 
uncertainties, we give a rough estimate of a mean-free 
path by assuming a monopole velocity v = a&v « 20 m/ s. 
Here, ad =4. 34 A denotes the distance of neighboring 
tetrahedra and v is the rate of single spin flips, which 
we roughly estimate by the monopole excitation energy 
E m » 2.2 K to v = E m /h « 4.6 • 10 10 s -1 . Assuming 
a temperature- independent v, i — 3D mag /v yields a lin- 
ear scaling from D mag (T) to £(T). As shown in figured] 
(right scale), within this estimate £(T) reaches the /jm 
range, i.e. ~ 1000 unit cells, for T < 500 mK, which may 
be understood from the low monopole density. However, 
even around ~ 2 K, I is still of the order ~ 100 unit cells 
although almost every second tetrahedron is in an ex- 
cited 3in-lout or lin-3out configuration. This suggests 
that at least towards higher temperatures, more complex 
hopping models have to be involved to describe K mag (T). 

In conclusion, we observe clear evidence for a consid- 
erable magnetic contribution n mag to the heat transport 
in the spin- ice material Dy 2 Ti 2 07. At constant tem- 
perature, the magnetic-field dependent decrease of k(B) 
correlates with the increase of the magnetization M (B), 
which measures the population of the particular 2in-2out 
configuration becoming the field-induced non-degenerate 
ground state. This reveals that the complete suppression 
of Kmag in a magnetic field of about 0.5 T arises from 
the loss of monopole mobility or, vice versa, the large 
«mag(0 T) is a consequence of the zero- field ground state 
degeneracy. Our data also reveal a strong increase of the 
relaxation times below about 600 mK. Including specific 
heat data, we derive a strong increase of the diffusion co- 
efficient below 1 K, which is most probably related to the 
fact that the monopole gas is highly dilute towards low 
temperature. In order to derive quantitative information 
about, e.g., the mean-free path or the velocities, theoret- 
ical models about the monopole dynamics are required. 

We acknowledge fruitful discussions with L. Fritz, 
M. Griininger, K. Kiefer, and B. Klemke and finan- 
cial support by the Deutsche Forschungsgemeinschaft via 
SFB 608. 



* E-mail: tl@ph2.uni-koeln.de 
[1] C. Castelnovo, R. Moessner, and S. L. Sondhi, Nature 

451, 42 (2008) 
[2] D. J. P. Morris et al, Science 326, 411 (2009) 
[3] S. R. Giblin et al, Nat. Phys. 7, 252 (2011) 
[4] C. Castelnovo, R. Moessner, and S. L. Sondhi, Phys. Rev. 

B 84, 144435 (2011) 
[5] D. Slobinsky et al, Phys. Rev. Lett. 105, 267205 (2010) 
[6] H. Kadowaki et al, J. Phys. Soc. Jpn. 78, 103706 (2009) 
[7] L. D. C. Jaubert and P. C. W. Holdsworth, J. Phys.: 

Condens. Matter 23, 164222 (2011) 
[8] S. T. Bramwell et al., Nature 461, 956 (2009) 
[9] S. J. Blundell, Phys. Rev. Lett. 108, 147601 (2012) 
[10] L. Yaraskavitch et al, Phys. Rev. B 85, 020410 (2012) 



■5 



[11] J. F. Nagle, J. Math. Phys. 7, 1484 (1966) 
[12] A. P. Ramirez et ai, Nature 399, 333 (1999) 
[13] S. T. Bramwell and M. J. Gingras, Science 294, 1495 
(2001) 

[14] Z. Hiroi et al, J. Phys. Soc. Jpn. 72, 411 (2003) 
[15] T. Sakakibara et al, Phys. Rev. Lett. 90, 207205 (2003) 
[16] K. Matsuhira et al, J. Phys.: Condens. Matter 14, L559 
(2002) 

[17] B. Klemke et al, J. Low Temp. Phys. 163, 345 (2011) 
[18] These differences hardly change the entropy j c/TdT. 
[19] S. Dunsiger et ai, Phys. Rev. Lett. 107, 207207 (2011) 
[20] Due to strongly increasing time scales below 500 mK 



and in order to avoid systematic errors induced by time- 
dependent experimental conditions, the heating curves 
Tp(t) have to be cut off after suitable times, typically 
5-10 minutes. Further relaxation can be estimated to 
increase c p by only a few percent. Thus, the presented 
data represent lower bounds for c p . 
[21] A. Anderson et al, Phys. Rev. 183, 546 (1969) 
[22] S. Saha et al, Phys. Rev. B 78, 214102 (2008) 
[23] N. Kamaraju et al, Phys. Rev. B 83, 134104 (2011) 
[24] In order to capture n in n(Q T, 400 mK, t) weaker stabi- 
lization criteria were used. 
[25] S. Erfanifam et al, Phys. Rev. B 84, 220404 (2011) 
[26] H. Aoki et al, J. Phys. Soc. Jpn. 73, 2851 (2004) 



